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Controlling entanglement sudden death in cavity QED by classical driving fields 
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We investigate the entanglement dynamics of a quantum system consisting of two-level atoms 
interacting with vacuum or thermal fields with classical driving fields. We find that the entanglement 
of the system can be improved by adjusting the classical driving field. The influence of the classical 
field and the purity of the initial state on the entanglement sudden death is also studied. It is 
shown that the time of entanglement sudden death can be controlled by the classical driving fields. 
Particularly, the entanglement sudden death phenomenon will disappear if the classical driving fields 
are strong enough. 
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I. INTRODUCTION 

Entanglement, one of the most striking features of 
quantum mechanics, has been considered as a key re- 
source of quantum information processing [J 0, S lj]- 
In recent years, the manipulation of quantum entangle- 
ment for the system of the cavity quantum electrody- 
namics(QEp) has been extensively investigated^, [a, 0, 
i, i, [IS 111 [ll [S [II- The cavity QED, where atoms 
interacting with quantized electromagnetic fields inside a 
cavity, can be used to create the entanglement between 
atoms in cavities and establish quantum communications 
between different optical cavities. It has been shown that 
the entanglement can arise in the interaction of a atom 
with a cavity field in a thermal state |6, 7] and two atoms 
with a cavity field [l3|, [IJ] . In Ref. [l5| , Solano et al. have 
shown that multipartite entanglement can be generated 
by putting several two-level atoms in a cavity of high 
quality factor. With the help of a strong classical driv- 
ing fields, the Schrodinger cat state and other entangled 
states can be produced. 

On the other hand, many efforts have been devoted to 
the study of the evolution of the joint system formed by 
two qubits [SE3, [H, [H, [H dU, [ll [H . Particularly, Yu 
and Eberly 16] pointed out that the single-qubit dynam- 
ics and the global dynamics of an entangled two-qubit 
system subjected to independent environments may be 
rather different. For example, for a single-qubit system 
subjected to an environment, the local coherence decays 
asymptotically. However, the entanglement of an entan- 
gled two-qubit system may disappear within a finite time 
during the dynamics evolution. The nonsmooth finite- 
time disappearance of entanglement is called "entangle- 
ment sudden death" (ESD). Recently, it is reported that 
the ESD phenomenon has been observed in a quantum 



optics experiment [2J]. 

In the present paper, we propose a scheme to enhance 
the entanglement of a quantum system consisting of two- 
level atoms within cavities by applying and controlling 
classical driving fields. The influence of the classical field 
and the purity of the initial state on the ESD is also stud- 
ied. First, we study the entanglement of the two-level 
atom and the field by employing logarithmic negativity 
[25l . l26l | . Then, we consider a quantum system consisting 
of two nonintcracting atoms each locally interacting with 
its own vacuum field. The two atoms, which are initially 
prepared in extended Werner-like (EWL) states [23| , are 
driven by two independent classical fields additionally. It 
is shown that the ESD may appear in this system and the 
ESD time can be retarded by increasing the purity of the 
initial state of the two atoms. In addition, the amount of 
the entanglement of the two atoms can be significantly 
increased by applying classical fields. We find that the 
time of entanglement sudden death can be controlled by 
the classical driving fields. It is interesting to point out 
that the entanglement sudden death phenomenon disap- 
pear if the classical driving fields are strong enough. 

The present paper is organized as follows. In section 
II, we propose a scheme to improve the entanglement 
of a two-level atom interacting with vacuum or thermal 
field by applying a classical driving field. In section III, 
we study the infiuence of the classical driving fields on 
the ESD when the atoms are initially prepared in EWL 
states. A conclusion is given in section IV. 



II. ENTANGLEMENT DYNAMICS OF A 

TWO-LEVEL ATOM IN A CAVITY WITH A 

CLASSICAL DRIVING FIELD 



The model 
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Now, we consider a system consisting of a two-level 
atom inside a single mode cavity. The atom is driven 



by a classical field additionally. The Hamiltonian of the 
system can be described by [15] 



+A[e " Uj^ + e cr_j, 



(1) 



where w, wo and Wc are the frequency of the cavity, atom 
and classical field, respectively. The operators a^ and a± 
are defined by cTz = \e){e\-\g){g\, a+ = \e){g\, and a. = 
cr\ where |e) and \g) are the excited and ground states 
of the atom. Here, a and a^ are the annihilation and 
creation operators of the cavity; g and A arc the coupling 
constants of the interactions of the atom with the cavity 
and with the classical driving field, respectively. Note 
that we have set h = 1 throughout this paper. 

In the rotating reference frame the Hamiltonian of the 
system is transformed to the Hamiltonian Hi under a 
unitary transformation Ui — exp (—iujctaz/'^) 



Hi = 



ulHUi-iUl^ 



r(i) 



= H\'' + H{ 



(2) 



with 



H 



(1) 



uoa) a + g(e*'^''*cr+a + e 



■'■'^■''a-a) 



H 



(2) 



Ai 



a^ + A(o-+ +0--), 



(2) 



(3) 



and Ai — loq — lJc- Using the method similar to that used 
in Ref. [27| , diagonalizing the Hamiltonian H\^ ' , and ne- 
glecting the terms which do not conserve energies (rotat- 
ing wave approximation), we can recast the Hamiltonian 
Hi as follows 



Hi 



ioa' a + 



J7i sin( 



0-+ +(T-] 



.9C0S^^[e'-* 



, sm y 2^ ■ 1^ \ , 1 /4\ 

x( 0-2 + cos -CT+ — sm -cr_)a + h.c\, (4) 



with 9 = arctan(-^). Here h.c stands for Hermi- 
tian conjugation. The Hamiltonian (4) can be diago- 
nalizcd by a final unitary transformation U2 with U2 = 
exp [^^(cr_|- -f- (T-)]. Then, we can rewrite the Hamilto- 
nian of the system 

t uj' sin 9 , sinfl 
H2 = u;a^a-\ {<t+ + a^) + g [{ —a^ 



-cos — 0"+ — sin —a-)a + h.c], 



(5) 



where w' = a/A^ + AX'^+Uc and g' = g cos^ | . It is worth 
noting that the unitary transformations Ui and U2 are 
both local unitary transformations. As we known the en- 
tanglement of a quantum system does not change under 
local unitary transformations [2^. Thus, the entangle- 
ment of the system considered here will not be changed 
by applying the local unitary transformations Ui and U2 ■ 



B. Entanglement dynamics 

We study the entanglement of the system by employ- 
ing the logarithmic negativity. For a bipartite system 
described by density matrix p, the logarithmic negativ- 
ity is 



Eip) = log2{l + 2N)^log,\\p^\\, 



(6) 



where p"^ is the partial transpose of p , \\p'^\\ stands for 
the trace norm of p'^ and N is negativity which is defined 
byJH 



N: 



\P 



1 



(7) 



which is just the absolute value of the sum of the nega- 
tive eigenvalues of p^. The logarithmic negativity has 
been proven to be an operational good entanglement 
measure [26|. We first assume that the atom is initially 
prepared in the state |-|-) = cos||e) -I- sin||5) and the 
field is in the Fock state \n). Then, we can find that the 
state vector at time t is 

\^{t)) = a„(i)|+,n)+/3„+i(t)|-,n+l), 
an{t) = e'^^*/2[cos(a.t)-^sin(f}„i)], 



Pnit) 



-ig cos^ -y/ne^''^^^/^ sin {nn-it)/nn-i, 



A2 = \/(wO - Wc)2 + 4A2 +UJC-UJ, 
n^ = yM + („+l)(5COs2^)2, 

1+) = cos-|e)-|-sin-|5), 

f) f) 

|-> = -sin-|e)-hcos-|5), 



(8) 



Inserting the state vector into Eq.(7) leads to the loga- 
rithmic negativity of the state |V'(i)) 



i?(|^»=log2[l + 2|a„(i)/3„+i(t)|]. 



(9) 



Next, we consider the entanglement dynamics of the 
system when the atom and the field are initially prepared 
in the state |-|-) and the thermal state, respectively. The 
initial density matrix of the system is 



p(0) = |+)(+|®(^P„|n)(n|), 



n=0 



m 



Pn 



(l + m)"+i' 



(10) 



where rn = l/{e^'^ — 1) is the mean photon number at 
the inverse temperature f3 and uj is the frequency of the 
optical mode. A straightforward calculation shows that 




FIG. 1: The logarithmic negativity E{\tp)) is plotted as a 
function of t with uj = bjiJo — 1 for cjc = A = O(dotted 
line) and uic = X = 2(solid line). Upper panel: The field is 
initially prepared in the vacuum state. Lower panel: The field 
is initially prepared in the Fock state |1). 




FIG. 2: The logarithmic negativity E{p) is plotted as a func- 
tion of t with u! — 5, too — 1 for lJc = A = O(dotted line) and 
ujc ~ X = 2(solid line). Upper panel: The mean photon num- 
ber of the thermal field is 0.1, i.e., m = 0.1. Lower panel: The 
mean photon number of the thermal field is 0.3, i.e., m = 0.3. 



the density matrix at time t is 



III. CONTROLLING ENTANGLEMENT 
SUDDEN DEATH 



Pit) - £{p„-i|/3n(i)n->(-|®|n)(n|+p„|a„(t)|^ 

n=Q 

x|+)(+| (E) \n){n\ + b„a„(t)/3*+i(i) 
x|+)(-|(g) \n){n+l\ + h.c]} 



fll) 



which leads to the logarithmic negativity of the above 
density matrix 

oo 

E{p) = log2[l + ^(|C|-C)], 

n=0 
C = 2^Pn-l\(3nf + P7i+l\an+lf 



-[(p„_i|/3„p -p„+i|q;„+iP)^ 
+4|p„a„/3:+ip]V2}. 



(12) 



The logarithmic negativity is plotted as a function of time 
t in Fig.l that when the field is initially prepared in the 
vacuum state |0) or in the Fock state |1). For the sake of 
simplicity, we set g — 1 in the numerical calculation. We 
can see clearly from Fig.l, the atom- field entanglement 
can be increased by applying the classical driving field. 

In Fig. 2, the logarithmic negativity E{p) is plotted as 
a function of t with uj = b,LUQ = I for Uc = X = (dotted 
line) and lUc = X = 2 (solid line). The mean photon num- 
ber m for the upper panel and the lower panel are 0.1 
and 0.3, respectively. It is not difficult to find that the 
entanglement of the two-atom system decreases with the 
mean photon number m. In other words, the amount of 
entanglement decreases with the increase of the temper- 
ature of the thermal field. However, we can increase the 
amount of entanglement by controlling the classical field. 



In the previous section, we have investigated the en- 
tanglement dynamics of a two- level atom interacting with 
a cavity field. Recently, Bellomo et al. [23, [^ investi- 
gated the influence of the Markovian and Non-Markovian 
effects on the dynamics of two-qubit entanglement. In 
this section, we consider a quantum system consisting of 
two noninteracting atoms each locally interacting with 
its own vacuum field. The two atoms, which are initially 
prepared in EWL states, are driven by two independent 
classical fields additionally. Our purpose here is to inves- 
tigate the influence of the classical driving field and the 
purity of the initial states on the entanglement dynamics 
of the system. 



A. Reduced density matrix of the two-atom system 

Here, we study the entanglement dynamics of two in- 
dependent atoms each locally interacting with a vacuum 
field by using the procedure of Ref. ^23| . In addition, each 
atom is driven by a classical field. In the following, we 
find the reduced density matrix of the two-atom system. 
In the basis |±), the initial state of the atom interacting 
with the vacuum field can be cast in the following form 



-(0) P, 
-(0) P- 



(0) 
(0) 



®|0)(0|, 



(13) 



where |-f) and |— ) are defined by Eq.(8). Using the state 
vector in Eq.(8) and tracing over the freedom of the field, 
we obtain the single-atom density matrix evolution 



P++W p+-{t) 

p_+(t) p__(i) 



(14) 



where 

p++{t) = |ao(i)lV-(0), 

p__(t) = [l-\ao(mp++{0), 

p+^{t) = ao{t)p+_{0) = p*M). 



(15) 



By making use of the above equations, we can find the 
dynamics of the two atoms via a purely algebraic way. 
Note that this procedure is applicable to arbitrary initial 
states of the whole system. In the basis |1) = |++), |2) — 
I + -), |3) = I - +)', |4) = I - -), and using Eq.(8) and 
Eq.(15), the density matrix for the two-atom system is 
calculated as follows 



Piiit) 

P22(t) 

Pssit) 

Pii{t) 

Pl2it) 
Pl3{t) 
Pl4{t) 
P23{t) 
P24{t) 
P34{t) 



|ao(i)lVii(0), 



|aoWl']Pii(0) 
\ao{t)\^]pii{0) 



\ao{t)\^P22{0), 
|ao(i)|'p33(0), 



|ao(i)Hl- 
\ao{t)f[l- 
[l-|ao(t)|']Vii(0) + P44(0) 

+ [l-|ao(t)P][p22(0)+P33(0)], 

ao{t)\ao{t)\^pi2{0), 

aoit)\aoit)\^pi3i0), 

aa{t)'^pu{0), 

I ao(i) 1^23(0), 

ao(t)[l - |ao(t)|']/5i3(0) + ao{t)p24{0), 

ao{t)[l - |ao(t)fl/9i2(0) + ao(i)/034(0), (16) 



with pij{t) = p*i{t). We would like to point out that the 
above procedure allows us to obtain the reduced density 
matrix of the two-atom system for any initial state. 



B. Extended Werner-like states and logarithmic 
negativity 

We assume the initial states of the two-atom system 
are the extended Werner-like states defined by 



p$(0) - r|<&)(<&| 



1-r 



4 ' 
p*(0) = r|vl/)(*| + i-^/, 






P\ 
P\ 






(17) 



where r is a real number which indicates the purity of 
initial states, / is a 4 x 4 identity matrix, p and v are 
complex numbers with |/zp-|-|j^p = l. It is noted that 
the EWL states belong to the class of the 'X' states. 
Explicitly, if the density matrix of a quantum state is of 
the form 



(18) 



Pll 








Pl4 





P22 


P23 








P*23 


P23 





Pl4 








P44 



The EWL states have the following advantages. First, 
we can easily find that the density matrix at arbitrary 
time t is still X structure under the single atom evolution 
determined by the Hamiltonian H2 (in the basis |1) = 
l + +),|2) ='! + ->, |3) = |-+),|4) = I--)) if the 
initial state is X structure, such as the Bell states and 
the EWL states. Second, the EWL states allow us to 
clearly show the influence of the purity and the amount 
of entanglement of the initial states on the entanglement 
dynamics simultaneously. The purity of the EWL states 
are dependent on the parameter r and the amount of the 
entanglement of the EWL states are related to r and p. 
If r = 1 the EWL states reduce to the Bell- like states |<i>) 
and I*). In the case of r = 1, ^ = z/ = l/\/2 the EWL 
states become the Bell states while in the case of r = 
they are the maximally mixed states. 

Note that we assume that the initial states of the sys- 
tem are the EWL states. The X structure of the EWL 
states is maintained during the evolution of the two-atom 
system. The explicit analytical expression of the logarith- 
mic negativity for the two-atom system can be obtained 
as 



Ex 

N 



log2(l + 2iV), 



max{0, -[\/{p22 - P33Y + 4|pi4p - P22 - P33]} 



- max{0, -Wipii - PAif +4|p23p 



-Pll - P44\]- 



(19) 



then it belongs to the class of the X states. 



We plot the logarithmic negativity as a function of time 
for different values of the purity of the initial states r and 
the strength of the classical field A. In what follows, we 
denote the entanglement of the two-atom system which 
is initially prepared in p$(0) (p.i((0)) by £■$ {Exy). 

The logarithmic negativity _£$ and E<i, arc plotted as a 
function of A and ujc with a; = o^o = 1, |mP = 1/6, \v\'^ — 
5/6, r = 2/3, t = 2 in the upper panel and the lower 
panel of Fig. 3, respectively. From Fig. 3, one can clearly 
see that if the parameters uj^ and A are small, then the 
logarithmic negativity E^ and E^^, are zero at time t = 2. 
The situation is different if the parameters of the classi- 
cal driving fields increase. When the values of Wc and A 
are lager enough the two atoms become entangled. The 
behavior of £^$ and E^i are different when the classical 
fields are applied. For example, the area for _E$ = is 
small than that of iJiji = which means that the entan- 
glement of the two-atom system is more sensitive with 
the classical driving fields if they are initially prepared 
in p$. This feature can be see more clearly in Fig. 4 and 
Fig. 5. It also interesting to note that the entanglement 
of the two-atom system can be significantly increased by 
controlling the classical driving fields both for the initial 
states p$ and pij,. 

In Fig. 4, we plot the logarithmic negativity of the two- 
atom system as a function of time t for several values 
of ijjc and A. The upper panel of Fig. 4 shows that the 
ESD occurs without the classical field (see the dotted 




FIG. 3: Upper panel: The logarithmic negativity _E$ is 
plotted as a function of A and Uc with ui = luq = l,|/ip ~ 
l/&,\v\^ = 5/6, r = 2/3, i = 2. Lower panel: The logarith- 
mic negativity £* is plotted as a function of A and ujc with 



■ Wo 



W 



1/6, 



= 5/6,r = 2/3,i = 2. 



line in this panel). After finite dark periods, the en- 
tanglement E^ revivals completely. The entanglement 
can be increased by applying the classical field(see the 
dashed and solid lines of the panel). It is interesting 
to point out that the ESD disappears when the classi- 
cal field is strong enough. Comparing the upper panel 
with the lower panel of Fig. 4, we find that the behavior 
of entanglement E<i, and E^, are different since the ESD 
appears in E^^j with ljc = 0.5, A = 1 while there is no ESD 
for E^ with the same values of lOc and A. Besides, the 
entanglement dark periods of E^, is longer than that of 
Ei^ when all of the other correspongding parameters are 
the same. In other words, it takes more time for E<^ to 
revive the initial entanglement. The parameters /i and v 
in Fig. 4 are chosen to be \/TJ& and -y/S/S, respectively. 
It is worth noting that different choices of /z and v do not 
give entanglement dynamics qualitatively different from 
the case treated here. 

Another aspect of interest is how the entanglement dy- 
namics is influenced by the purity of the initial states. 
In order to show it intuitively, we plot the logarithmic 



0.4 -■ 
0.3- 
lil* 0.2- 
0.1- 
0.0 — 




FIG. 4: Upper panel: The logarithmic negativity i5<j is plot- 
ted as a function oft with |/i| = 1/6, \v\ = 5/6, r = 2/3, to = 
u)o = 1 for LJc = A = (dotted line), Uc = 0.5, A = 1 
(dashed line), and uOc = 0.5, A = 2 (solid line). Lower panel: 
The logarithmic negativity E^a is plotted as a function of 
t with |/.ip = 1/6, |z/p = 5/6, r- = 2/3, c^ = o^o = 1 for 
ctJc = A = (dotted line), ujc = 0.5, A = 1 (dashed line), 
and Uc = 0.5, A = 2 (solid line). 





FIG. 5: Upper panel: The logarithmic negativity i5<j is plot- 
ted as a function of t with |/ip = 1/6,11^1^ — 5/6, a; = 
ujo = l,LOc = 0.5, A — 1 for r = 0.5 (dotted line), r = 0.7 
(dashed line), and r = 1 (solid line). Lower panel: The 
logarithmic negativity Eif is plotted as a function of t with 
I^P = 1/6, |i/p = 5/6, Lj^LJo = l,iOc = 0.5, A = 1 for r = 0.5 
(dotted line), r = 0.7 (dashed line), and r — 1 (solid line). 



negativity of the two-atom system as a function of t for 
different values of r in Fig. 5, fixing the parameters jj. and 
i/. From Fig. 5, one can see that the entanglement dark 
periods increase with the decrease of the purity of the ini- 
tial states. The classical driving field again can increase 
the amount of the entanglement of the system. 



IV. CONCLUSIONS 

In the present paper, we have considered a quantum 
system consisting of one two-level atom interacting with a 
single mode field. The atom is driven by a classical field 
additionally. The entanglement dynamics is described 
by the logarithmic negativity. We first investigated the 
entanglement dynamics of the atom-field system. The 
field is initially prepared in the Fock state or the thermal 
state. We find that the atom-field entanglement can be 
increased by applying the classical driving field. 

Then, we consider a quantum system consisting of 
two noninteracting atoms each locally interacting with 
its own vacuum field. The two atoms, which are driven 
by two independent classical fields, are initially prepared 
in EWL states, i.e., p$ and p$. We find that classical 
driving fields can increase the amount of entanglement of 
the two-atom system. The behavior of entanglement i?$ 



and E^i, are different when the classical fields are applied. 
The entanglement of the two-atom system is more sen- 
sitive with the classical driving fields if they are initially 
prepared in p$ . It is worth noting that the time of entan- 
glement sudden death can be controlled by the classical 
driving fields. Specially, the entanglement sudden death 
phenomenon will disappear if the classical driving fields 
are strong enough. It is interesting to verify the scheme 
of our paper in cavity QED experimentally. 
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